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Abstract 

Inspired by the corresponding problem in QCD, we determine the pressure of massless O(A^) 
scalar field theory up to order in the weak-coupling expansion, where g^ denotes the quartic 
coupling constant. This necessitates the computation of all 4-loop vacuum graphs at a finite 
temperature: by making use of methods developed by Arnold and Zhai at 3-loop level, we 
demonstrate that this task is manageable at least if one restricts to computing the logarithmic 
terms analytically, while handling the "constant" 4-loop contributions numerically. We also 
inspect the numerical convergence of the weak-coupling expansion after the inclusion of the 
new terms. Finally, we point out that while the present computation introduces strategies 
that should be helpful for the full 4-loop computation on the QCD-side, it also highlights 
the need to develop novel computational techniques, in order to be able to complete this 
formidable task in a systematic fashion. 
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1. Introduction 



Motivated for instance by hydrodynamic studies of heavy ion collision experiments, and dark 
matter relic density computations in cosmology, a lot of theoretical work has been devoted 
to the perturbative determination of the pressure of hot QCD in recent years. As a result 
of 3-loop and 4-loop computations, corrections to the non- interacting Stefan-Boltzmann law 
have been determined up to relative orders 0{g'^) 0{g^) [21 E]' C){g^ln{l/g)) [5], 
where g denotes the renormalized strong coupling constant. The first presently unknown 
order, 0{g^), contains a non-perturbative coefficient ^l6j, but that can also be estimated 
numerically O |8] . All orders of g are available in the formal limit of large iVf [9] , where A'^f 
counts the number of massless quark flavours. Similar results have also been obtained for the 
case of non-zero quark chemical potentials at finite [10] and large [11] Nf. Moreover, first steps 
towards the inclusion of finite quark masses, important for phenomenological applications, 
have been taken [12]. Finally, coefficients up to the order 0{g^) are available even for the 
standard electroweak theory, at temperatures higher than the electroweak scale [13] . 

Conceptually, it would be quite desirable to extend these results up to the full order 0{g^). 
The reason is that this is the first order where at least the leading contributions from the 
various momentum scales relevant for hot QCD, k ~ 2ttT, gT, g'^T, have been fully accounted 
for. As mentioned, the non-perturbative input needed to describe the effect of the softest 
momenta k ~ g^T is already available [71[8]. Moreover, several perturbative contributions of 
0{g^) are known: in the notation of Ref. [1], /3m |14j (which accounts for the scales k ~ gT), 
(3e2 [E], Pe3 [E], as well as /3e4 and /3e5 [IZ] have been computed. Nevertheless, a single 
coefficient, /3ei, coming from the hard scales k ~ 2ttT, remains undetermined. 

The recipe to compute Pei is in principle simple: it is defined to be the "naive" (i.e. 
unresummed) 4-loop contribution to the pressure, computed by regulating all divergences 
via dimensional regularization in d = 4 — 2e dimensions. This computation contains infra- 
red (IR) divergences, which manifest themselves as an uncancelled 1/e-divergence in the 
final (renormalized) result. This IR pole gets only cancelled once the contributions of the 
soft modes are properly resummed, a step that has already been completed, and produces 
another 1/e-pole [1], with the opposite sign. 

Unfortunately, the practical implementation of this computation is far from trivial. In fact, 
to show that it is feasible at all, it is the purpose of the present paper to demonstrate that 
the 4-loop order can at least be reached in models somewhat simpler than QCD. 

The model that we consider is scalar field theory, with a global 0(A) symmetry, i.e. the 
"A(/)^" -theory. To keep the analogy with QCD in mind, we follow here the frequent convention 
of denoting X = g^. The pressure of this theory has been computed to high orders in g in 
parallel with that for QCD: the orders 0{g^) [lEld], Oig^) [El [20], and 0{g^ln{l/g)) [20] 
have been reached. The order 0{g^) involves a coefficient analogous to /3ei, which again 
contains an uncancelled 1/e-pole. It is the goal of the present paper to compute this /?ei (as 
well as all other relevant coefficients), showing that the pole cancels, and thus to determine 
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the full pressure up to 0{g^), at any finite N . 

It is perhaps worth stressing that even though the technical challenges addressed in this 
paper have a direct counterpart in QCD, there is of course the conceptual difference that in 
our model the order 0{g^) contains no non-perturbative coefficients, since the scale k ~ g^T 
does not exist in scalar field theory. Therefore no lattice studies of the type in Refs. [i, ^ 
need to be invoked here. 

The plan of this paper is the following. We start by carrying out the naive 4-loop computa- 
tion of the pressure in Sec. [21 leading to the determination of the coefficient /3ei • In Sec. O we 
elaborate on how the naive computation can be re-interpreted and incorporated in a proper 
setting such that all divergences cancel. This leads to our final finite result. We discuss 
various formal and numerical aspects of the result in Sec. IH and conclude in Sec. [5l 



2. Naive 4-loop computation 

Our starting point is the bare theory 

4-2e N 1^2 

= 2 E E ^,<l^^9^^<t>i + j/b^^' (E 't'^^) ' (2.1) 

^1=1 i=l 1=1 

where A is the scale parameter introduced in connection with dimensional regularization. We 
work in Euclidean metric throughout. The fields (pi have the dimension [GeV]^~^, while the 
bare coupling g^ is dimensionless. The theory can be renormalized in the MS scheme by 
setting 

«l=»' + 7i4^ + 7i^f4 + (2.2) 



(4^)2 e (47r)4Ve2 2e 

where (see, e.g., Ref. [21]) 



D D 

The observable that we are interested in, is the pressure, or minus the free energy density, 
of this theory, as function of the temperature T and the renormalized coupling constant g. 
Formally, the pressure is given by the path integral 

p{T) = hm ^ In / VcPi exp(-5B) , (2.4) 

V^oo V J 

where V is the volume, V = /d'^~^'^x, and is the bare action, Sb = Jq^^t Jd^~^'^xCB ■ 
The temperature T is given by T = 1/(3, and the path integral is taken with the usual 
periodic boundary conditions over the r-direction. Note that in the presence of dimensional 
regularization, the pressure thus defined has the dimension [GeV]'^"^'^. 
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Figure 1: T/ie i-, 2-, 3- and 4-loop full theory diagrams contributing to the pressure. 
Going into momentum space, we define the usual (sum-)integrals as 



(2vr)3-2e \ AtT JJ (2^)3-2^ 

Tj;/, (2.6) 



PO 

where po stands for bosonic Matsubara momenta, po = 27rnT, n E Z. Furthermore, A is the 
MS scale parameter, while A is the MS one; these are related through A^ = 47rA^ exp(— 7e). 
We then define 

^ iow^ ■ <^-'> 

Up to 4-loop order, the strict loop expansion for the pressure of our theory contains the 
graphs (a)-(h) in Fig. [H which we denote by the symbols larlh- Denoting furthermore 
Xn = 2^, and employing the sum-integrals defined in Eqs. (|2.7p - (j2.9p . we can immediately 
write down their expressions in the forms: 

^2 



4 = 


Ar^r4[l + 0(6)], 


(2.10) 


h = 


iV(iV + 2) 2,^,2 

24 5b (2^1 ) > 


(2.11) 


Ic = 




(2.12) 


Id = 




(2.13) 


h = 




(2.14) 


h = 


N{N + 2)3 

1296 ' 


(2.15) 
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N{N + 2){N + 8) g 



1296 



(2.16) 
(2.17) 



Taking now into account the renormalization of the couphng constant up to 2- loop order, we 
must replace the bare coupling by the renormalized one, g'^, through Eq. (|2.2p . Recalling 
that a weak-coupling expansion (i.e., an expansion in g'^) of the path-integral in Eq. (j2.4p 
does not coincide with the loop expansion at finite temperatures, because of well-known IR 
divergences in multiloop graphs, we will denote the result for the sum of these graphs in the 
following by pe{T), in contrast to the physical pressure denoted by p{T); for the latter we 
assume that a consistent evaluation has been carried out to a certain order in g'^, irrespective 
of how many loop orders this takes. Summing together the graphs shown leads then to the 
following expression for pE (T) : 



N-T^[l + 0{e)]-g'x 



24 



(N + 2)l2 



N + 8 
(4^ 



+ 



1296 

2{N + 2){N + 8) 



{N + 2)^Iifl3 + ^{Ii'^ 



I2 + 



1 



3 
2 

{N + 8f 3(3iV + 14) 



{N + 2)2 (T2) 



(47r)2e ' (47r)4 

+6(A^ + 2)li4 n(P)n(P) + {N + 



[HP)] 



3 3(iV + 



(47r)2e 



mp)x 



(2.18) 



Prom Ref. [T], we can immediately read off the results for most of the sum-integrals above, 
utilizing the formulae 



■^m 2'"~2™+l^m~2n+3/2 jnm— 2n+4 

2 1 



Yr(n-3/2 + e) 



vr 



J'2 



[ii{p)r 



(4vr) 



2^2 N 2 



12 



r(n) 
C'(-i) 



^+36 1n^+48^^_^ 



C(2n - m - 3 + 2e) , 

C(-3) 5 ^ ' 



(2.19) 



This leaves us with the task of evaluating the new 3- and 4-loop sum-integrals 



Si ^ 

^2 = 



n(P)n(p) , 



f^[mp)Y 



(47r)2e 



[n{p)Y 



(2.20) 

(2.21) 
(2.22) 
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where we have grouped together two terms in 5*2 for computational convenience. 

The sum- integrals in Eqs. (j2.2ip . (j2.22p can be evaluated through a very tedious if in 
principle straightforward application of the procedures and techniques that were pioneered 
by Arnold and Zhai in Ref. [T]. A detailed explanation of the steps that we have taken can 
be found in Appendix A. Here we simply quote the final results of that analysis: 



Si 



+ 



1 



1 



8(47r) 



In 



A^ 
47rT2 



+ 



3 In 
17 



A2 17 ^ C(-l) 
4^ + y + 7e + 



+ 



C'(-i) 



In 



A2 

+48.797635359976(4) \ +0{e) 



(2.23) 



rp4 



16(47r)4 
A2 ■ 



2 In 



A2 



1) 



10 X'{ 

4^+y-'^-+^cM) 



+ 



In 



47rT2 

rp4 



+ 



27E + 4 



C(-3) 
C(-3) 



In 



A2 



47rT2 



25.705543194(2) I 



512(47r)2 



+ 28.92504950930(1) } + 0{e) 



(2.24) 



In 52, we have on purpose separated one contribution (coming from the part denoted by S'g^ 
in Appendix A, cf. Eq. ()A.57P ) from the rest, as it contains an IR singularity, which will not 
be cancelled by the renormalization of the coupling constant in the full theory, but only by 
the ultraviolet (UV) singularities originating from the "soft" contributions to the pressure 
(cf. next section). The numbers in parentheses estimate the numerical uncertainties of the 
last digits shown. 

In order to present the full result, we introduce the following notation for the contributions 
of various orders to pE- 



OEi + g aE2 + 



9 



(4vr) 



9 



r(3Ei + 0{g\e) 



(2.25) 



Here the coefficients OEi, /3ei have been defined in analogy with Ref. [3]. They are dimension- 
less functions of the temperature and of the regularization scale. Inserting Eqs. (|2.19p . (j2.20p . 
(j2.23p and (j2.24p into Eq. ()2.18p . and using everywhere the MS scheme scale parameter A, 
we obtain the following expressions for these coefficients: 



"El 



aE2 



"ES 



90 ' 

N{N + 2) 

3456 
iV(iV + 2) 

10368 



{N + 8) In 



A 



1 



31 C'{ 
^ + (iV + 2),. + - + 12i^ 



■ C^(-3) 
'C(-3) 



(2.26) 
(2.27) 
(2.28) 
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/3ei 



N{N + 2){N + 8) vr^ 
41472 7 
N{N + 2 



31104 



+ 16iV 



+ 64 



107 



{N + t 
57e 



In 



80 
353 



+ 



AY 

Attt) ^ 
2 C(- 



27Eln- 



A 



^^^^ 36 



4 



In- 



A 



160 



+ 



7e 



3vr2 C'(- 



-11 



AttT 

1) 4 c(-3) y 4^r 

3C'(-3)\. A 



+ 



3.9753932807(2) 



+ 



1) 2 C(-3) J 47rr 



In- 



10.6970470860(3) 



(2.29) 



The expressions for uei, aE2, ctES agree with the results of Ref. [20] for = 1. 

We note from the first hne of Eq. (|2.29|) that the fully renormalized result indeed contains 
an uncancelled 1/e pole. Therefore, something must be wrong with the naive computation 
that we have carried out. We now turn to the correct procedure for determining p{T) up to 



3. Resummed 4-loop computation 

The reason that the computation carried out in the previous section leads to a divergent result, 
is that it ignores the fact that certain subsets of higher loop graphs amount to generating 
radiatively a mass for the fields (pi. In the presence of a mass, the results for IR sensitive loop 
integrals change. This cures the IR problem and leads to the correct weak-coupling expansion. 
A systematic way to implement such a mass resummation (as well as other resummations, for 
instance for the quartic coupling) goes via effective field theory methods, as we now review. 

At high temperatures, the relevant effective field theory framework is that of dimensional 
reduction [221 [23] . The basic observation is that the non-zero Matsubara modes are heavy, 
and certainly cannot cause any IR divergences. Thereby the computation of all static ther- 
modynamic observables can be factorised into two parts: to the contribution from the "hard" 
momentum modes, k ~ 27rT, and from the "soft" modes, k ~ gT (the scale k ~ g^T does 
not appear in scalar field theory) . For the pressure we will denote the two parts as phard and 
Psoft ) respectively. The effective theory determining Ps^ft reads 

^^ = 2^^ 9,-0,9j(/>i + - m| ^i^i + -glA^' (j^ c^i^i^ + . . . . (3.1) 

j=l i=l i=l i=l 

A series of infinitely many higher dimensional operators has been truncated given that, as 
power-counting arguments show (cf. Ref. [4]), they cannot contribute to the pressure at 
0{g^). The theory in Eq. ()3.ip describes the dynamics of the Matsubara zero modes, and 
thus lives in three dimensions. The parameters here are to be understood as bare parameters. 
In dimensional regularization, the dimension of (pi is [GeV]^/^~^ and that of g| is [GeV]^. 
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Note that for simplicity we have used the same notation for the fields in Eqs. (j2.ip . (j3.ip . 
even though they are independent integration variables, with a different mass dimension. 
The use of the subscript E is meant to keep in mind the analogy with the effective theory 
called EQCD in the context of QCD [3]. The determination of the effective parameters (or 
"matching coefficients") in Eq. ()3.ip up to 2-loop level dates back to Ref. [24], where the 
symmetry-breaking phase transition in the case of a massive scalar field was considered; the 
application of this theory to the computation of the pressure at very high temperatures, 
where the zero-temperature scalar mass can be ignored, as is done in this paper, was first 
pursued systematically in Ref. [20j, as far as we know. 

The factorization statement now reads that the physical pressure can be written as [20] 



p{T) =p^{T) +pyi{T) =Phard(r) + Psoit{T) , 



(3.2) 



where 



Vm{T) 



lim ^In / V(j)i exp(-S'E) 

V^QO V J 



(3.3) 



and Se = J d^~^'^xi2E • We adopt a notation in the following whereby the matching coefficient 
Pe is a bare quantity, like m| and g^, while Phard is defined as its MS scheme version. Similarly, 
Psoft is defined to be the MS version of pM- 

Now, if we were not to carry out any resummation — that is, if the mass parameter m| 
in Eq. (|3.ip were ignored as was the case in the previous section — then the path integral in 
Eq. (j3.3p would vanish order by order in dimensional regularization, because the propagators 
appearing in the computation would contain no mass scales. Therefore, according to Eq. (j3.2p . 
the computation of the full pressure without any resummation, produces directly the function 
Pe- In other words, the proper interpretation for the result of the previous section is to treat 
Pe{T) in Eq. (|2.25p as an UV matching coefficient: it contains the contributions to the physical 
pressure from the hard modes, k ~ 27rT. Interpreted this way, it is IR finite, because the 
soft contribution pm{T) has been subtracted (which happens automatically in dimensional 
regularization by ignoring the mass parameter |20j ) . 

Our task in the remainder of this section is to properly compute pM, as defined by Eq. ()3.3p . 
This result can be extracted directly from Ref. [H]. Setting g ^ 0, X g^/6, and (Ia ^ N 
there, we obtain 



A2' 



PM 



m| N 

4^y 



+ 



ff|m| N{N + 2) 
(47r)2 24 

g^mE N{N + 2) 
(47r)3 144 
ffl iV(iV + 2) 
(47r)4 1728 



1 + 2 



+ 



iN + 2) 



1 A2 

- +41n ^ +4-41n2 

e Ami, 



{N + 2) 



7 



(i) 



(ii) 



(iii) 



e 



(iv) 



Figure 2: The 1-loop and 2-loop graphs needed for determining the matching coefficients m\ 
and g\ in Eqs. h3. 5\) . h'j.b]) . respectively. 



1 A2 a/o\ 

- + 41n— ^ + 2 + 41n2-84^ 



Arn^ 



+ 0{e)] +0(^) . (3.4) 
J VmE/ 



We next have to determine the values of m| and g| that appear in Eq. (|3.4p . Following 
the notation in Ref. [4], these matching coefficients can be written as 

^4 



5| 



aE4 + OEse + 0{e 

,4 



+ 



g 



(47r)2 



(47r)2 
aE7 + /?E3e + ©(e^ 



«E6 + /3E2e + 0(e')] , (3.5) 

0(5')} • (3.6) 



Up to the orders indicated they are produced by the graphs in Fig. [2l Using methods 
explained in some detail in Ref. [16], we obtain the values 
A^ + 2 



aE4 

"ES 
Q:E6 
«E7 
Pe2 



72 

A^ + 2 

36 
N + 2 



In + 1 + ^ ^ ' 



4ttT 



C(-i) 



72 

iV + 8 
3~ 
iV + 2 



1 iV/ A \ 4/ A 



3 7e 3 C'{-1) 

AttT 2 2 2 C(-l) 



\ A 

In — ^ + 7e 



AttT 



(3.7) 
(3.8) 
(3.9) 
(3.10) 



1728 t 



24iV In^ 



16N 



64 



1 + 27e + 
1-7E + 



A 
4^ 

c(-i; 

C'(-i) 



+ 



11 



A vr^ C' 
^"4^+16+^^ + ^^CFT 



71 



C(-i) 

3 C^(-l) 3 r(-l) 
2 C(-l) 4 C(-l) 



7e 



In ^ 3 vr 

^4^^ 2 ^ 16 2 2 C(-l 



7eC'(-1) 71 
2 



Pes 



N + 8 



In' 



A 
4^ 



+ 27e In 



A 
4^ 



+ 



vr 



271 



(3.11) 
(3.12) 



Here 71 is a Stieltjes constant, defined through the series C(s) = l/(s — l) + X]J^o ')'"(" ^)"(^~ 
l)"/n!. The values of aE4, OEe agree with the results of Ref. [20] for = 1. The value of oe? 
agrees with the results of Ref. [24J, for = 1, 2, 40 In the following, we denote by 0£^f,l3^^ 



""^There is a sign error in Eq. (60) of Ref. [24], relevant for the case A*' = 1. 
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couplings from which the 1/e-divergences have been subtracted; this is in fact relevant only 
for aE6 and /?ei, cf. Eqs. I^M), 

We note from Eqs. (j3.5p . (j3.9p that the mass parameter has a divergent part. Prom 
the point of view of the effective theory, this divergence acts as a counterterm, 



(5m| 



gij^A AT + 2 
(4^)2 72e 



(3.13) 



Indeed, this agrees with the counterterm that can be determined within the effective theory, 
by just requiring renormalizability of |24j . It will be convenient for the following to also 
define finite parameters from which the divergence as well as terms proportional to e have 
been subtracted; we denote these by 



mi(A) 



"^3 (A) 

rp2 



9l 



g aE4 + 



9 „ MS 



(47r)^ 



(4^)' 



(3.14) 
(3.15) 



The parameter is renormalization group (RG) invariant up to the order computed, while 
m|(A) has dependence at order 5^, because the counterterm in Eq. (j3.13p has been subtracted. 



Re-expanding now the bare expression of Eq. (j3.4p by treating the counterterm in Eq. (|3.13p 
as a perturbation, we obtain a "renormalized" expression for pM- 



where, after sending e ^ 0, 
m|(A) N 



56 Ar(Ar + 2)(iV + 8)7r2 A^^Psoft 



(47r) 



41472e 



+ 



2^4 



(3.16) 



Psoft 
2^4 



+ 



An _3 
g|mi(A)iV(iV + 2) 

(47r)2 24 
glms{A)N{N + 2) 
(47r)3 288 
5f iV(A^ + 2)| 



;ln- 



A 



m3(A) 



+ 10-161n2-A^ 



+ 



{Air 



(N + Sf- 



864 
2 



4 In 



(iV + 2) 
A 



21n- 



A 



+ 1 - 41n2 



m3(A) 



+ 1 - 2 In 2 - 42 



C(3) 



+ 



(iV + 2) 



_JL_ 



(3.17) 



We note from Eq. (13.161) that the only divergence appearing in pM after the re-expansion of 
m| is of 0{g^). Therefore, none of the coefficients Qes , /3e2 , /3e3 in Eqs. (13. 5p . (|3.6p . which 
multiply terms of order ©(e), play a role in psoft at 0{g^). In this respect, the present theory 
differs from QCD, where the corresponding coefficients do play a role; the reason for the 
difference is that in QCD there is an 1/e-divergence in the second term of Eq. (13. 4p . 
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To conclude this section, we remark that it may be convenient, fohowing Ref. [T2j, to also 
express pE in terms of the parameter ^| . This implements a certain (arbitrary) resummation; 
the practical effect is minor (in fact, for the scale choice that we will make in Eq. (I4.16p . which 
leads to the vanishing of oey, there is no effect at all, cf. Eq. (I3.19P ). and we present this 
last step only in order to allow for a more compact numerical handling of the result. From 
Eqs. ^nU\f . ^M), we obtain 



A'^'PE _ 9^ N{N + 2){N + 8)7T^ A'^d 
T4 (47r)4 41472e 

where, after setting e — > 0, the finite function Phard can be written as 



(3.18) 



Phard _ 
2^4 



■ 



^4 

"El + 53«E2 + («E3 - aE2aE7) + 

JL 

(Any 



13^^ + 2aE2a|7 - 2aE3aE7] + 0(5, 



(3.19) 



In the numerical results of the next section, we refer to the various orders of the weak- 
coupling expansion according to the power of rh^,g^ that appear, with the rule 0(7713) = 



0(^3) = 0{g). In other words, ''0{g'^f denotes 0{g. 



n—k^ 



in the expression constituted 



by the sum of (|3.17p . (13.19p . If rh^,g^ were to be re-expanded in terms of g, one would 
recover the strict weak-coupling expansion (given in Eq. ()4.ip below); however, it is useful to 
keep the result in an unexpanded form, because this makes it more manageable, and because 
the unexpanded form introduces resummations of higher order contributions which may be 
numerically significant for the slowly convergent part Psoft in Eq. (j3.17p [Ml [SOI 13 HI [25] (in 
practice, though, the effects caused by this resummation are not dramatic). 



4. Results and discussion 



Inserting Eqs. (jS.Sp . (13. 6p into Eq. (13. 4p . summing together with Eq. (I2.25p . and sending 
e — > 0, we obtain the strict weak-coupling expansion for the pressure of our theory: 



P{T) 

rp4 



+ 

+ 

+ 
+ 



"El + 9 aE2 + 
N 3/2 



Att 3 
9' 
(47r)2 



-a 



E4 



+ 



N{N + 2) 

Q^E3 Cl!E4 



24 



/ N 1/2 

«E4 



(47r)3 2 

9' 



"E6 



fOMS 

Pei 



A^ + 2 
144 
A^(Ar + 2) 



+ 

81n- 



A 



gT^/oEi 



+ 10- 161n2 - N 



+ 



24 



a^g -I- aE4aE7 
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N + 2 
36 



21n- 



A 



216 



864 



A 



+ 1 - 2 In 2 - 42 



+ 0{g' 



(4.1) 



We note that the 1/e-divergences in Eqs. (I3.16p . (|3.18p have cancelled against each other, 
as must be the case for a consistently computed physical quantity. Inserting the expansions 
from Eqs. ()2.26p - (l2.29p and Eqs. (I3.7p ~ ()3.10p . finally yields the explicit expression 
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A number of simple crosschecks can be made on Eqs. (I4.2p -( 1T9]) . Making use of the RG 
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equation, 
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it is easy to verify that the result is RG invariant up to the order computed. Setting = 1, 
terms up to 0{g^) agree with Refs. [U [HI [20]. Finally, taking the limit N ^ oo with g'^N 
fixed, we get 



pjT) 

rp4 



+ 



90 

_5_ 
36 

6( ^ 

108 



12 



In 



A 



+ 7e 



V 9 

6^ +g' 



+ 



ln-^1 +2(7E-6)ln 



9\/2 
A 



+ 



AttT 



+ 6-127^ + 71 + ^ 



(4.11) 



where Q = gyN /Att. This agrees with Eq. (5.8) of Ref. [27]; unfortunately most of the non- 
trivial structures, like logarithms of g or the genuine 4-loop sum- integrals that we were only 
able to determine numerically, disappear in the large- A/^ limit. 

Let us finally evaluate our result numerically. Though the effect is moderate in practice, 
we reiterate that we find it convenient not to use Eq. ()4.ip for the numerical evaluation, but 
the unexpanded expression phard +Psoftj defined as a sum of Eqs. (j3.17p . (jS.lOp . 

If we want to present the numerical results as a function of T, we first have to insert 
a value for the renormalized quartic coupling g^ as a function of the scale A, because the 
temperature-dependence emerges in connection with logarithms related to the running of g"^ 
(cf. Eqs. ([12])-([13])). Defining hi = |/3i|/(47r)2 and &2 = |/32|/(47^)^ where /3i,/?2 are from 
Eq. (12. 3p . as well as q = 62/^1 and t = 2bi ln(ALandau/A), the 2-loop RG-equation reads 



dg^ 
dt 



-g' + ag^ 



(4.12) 



This equation can be solved up to a boundary condition. In formal analogy with QCD, we 
define the boundary condition such that 
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The solution then reduces to the equation 
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which for small A yields the approximate behaviour 
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Figure 3: The effective gauge coupling in Eq. 113. 15\) . as function ofT/Aiandau, where A-iandau 
is defined in Eq. j^. From left to right, N = 1,4, 16. The scale A has been varied within 
the range (0.5...2.0)Aopj (the grey band), with Aopt defined in Eq. (JA^, but the effects are 
almost invisible on our resolution. 



Moreover, following Ref. [26], we define an "optimal" scale according to the simultaneous 
1-loop "fastest apparent convergence" and "principal of minimal sensitivity" point obtained 
for the effective coupling g^, defined in Eq. (j3.15p : 

Aopt = 4^e-^^T . (4.16) 

The scale A will then be varied in the range (0.5... 2.0) Aopt around this point. Note that the 
scale choice in Eq. (|4.16p leads also to a formal simplification of the pressure; for instance 
the expression in Eq. (j4.1ip obtains the form 

p{T) _ ^^^li ^ 
T^ ^ ~9(r\ ~ 12 

We now first plot as a function of T/ALandau and N . The results are shown in Fig. [3l 
The scale A is chosen according to Eq. ()4.16p but the dependence on the choice is so small 
within the range mentioned that it is almost invisible. 

Finally we plot the pressure, normalised to the free result, in Fig. HI We have used an 
overall resolution as would be relevant for QCD, where one hopes to reach an accuracy on 
the 10% level or so; the inserts display finer structures that are invisible on this resolution 
but may still be of some academic interest. Since T/ALandau is a quantity with which it is 
difficult to associate anything in QCD, we choose to use g^ as the horizontal axis in this 
figure. If desired, the conversion to T/ALandau goes through Fig. [3l 

It can be seen that the patterns in Fig. |4]are somewhat similar to those familiar from QCD: 
the order 0{g^) strongly overshoots the free result, and the subsequent orders then slowly 
converge towards a common value. However, for N = 16, i.e. the case with the same number 
of scalar degrees of freedom as QCD, control is certainly lost once g^ ^ 4.0. It is perhaps also 
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Figure 4: The resummed perturhative result, p(T) = Phard{T) + Psoft(T) , normalised to the free 
Stefan- Boltzmann result denoted by psB, cls a function of the effective gauge coupling shown 
in Fig. 0. From left to right, N = 1,4, 16. For the contribution of 0{g^), the scale A has 
been varied within the range (0.5...2.0)Aopi (the grey band), with Kopt defined in Eq. ^.16^ . 

worth stressing that even at = 16, the results deviate fairly substantially from Eq. ()4.17p : 
relative corrections to the large- A^ limit can be as large as ~ 8/A^, cf. Eq. (j4.9p . 

A perhaps more satisfying view on the result can be obtained if the hard and the soft 
contributions to the pressure, Phard and Psofti are plotted separately. This has been done in 
Fig. [5] for N = 16, where the absolute values of each new order are shown. It is clear that at 
least the expansion for Phard does appear to converge up to g^ ~ 5.0; for psoft the situation 
is worse, since the 0{g^) result is larger than the 0{g^) result even at fairly small g^. This 
is, however, due to the fact that the 0{g^) result crosses zero at g^ ~ 2.5, and does not 
necessarily signal a total breakdown of the series. 

5. Conclusions 

We have computed in this paper the pressure of O(A^) scalar field theory up to order g^ in the 
weak-coupling expansion. In terms of the loop expansion, this corresponds to the inclusion 
of all 4-loop diagrams, as well as infinite subsets of higher-loop diagrams needed in order to 
cancel the infrared divergence of the naive perturbative computation. 

The main motivation for this paper has been "technological" : we have demonstrated that 4- 
loop sum-integrals are doable, with divergent parts and logarithms handled analytically, and 
constant parts evaluated numerically. The essential ingredients allowing for the evaluation of 
the single genuine 4-loop sum-integral that appeared in our computation, denoted by 52, were 
the realization that renormalizing the theory before carrying out the sum-integral, allows to 
simplify the structure that needs to be considered (cf. Eq. (j2.22j) ). as well as an application 
of mixed coordinate and momentum-space techniques (cf. Appendix A. 3). 

With further work, the constant parts of the 4-loop pressure (the numerical values in 
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N= 16 



N= 16 




Figure 5: The absolute values of the various order contributions to p hard/ O^ft) o-'^d Psoft/T^ 
(right), as a function of the effective gauge coupling shown in Fig. 0, for N = 16. The 
convergence appears to be better for p hard/ T"^ , but psoft/T"^ shows some convergence as well. 



Eq. (|4.9p ) might also be doable analytically. This, however, is not necessary from the QCD 
point of view, where the 4-loop contribution in any case involves a non-perturbative term OEj 
that can only be determined numerically [71[8]. 

The complete result, up to 0{g^), is shown in Eqs. (j4.2p ~ (j4.9p . A more compact form, re- 
summing hard contributions, is given by a combination of Eqs. (j3.17p . (j3.19p . with coefficients 
given in Eqs. (|377D - (|3T0]) and (|2:26D - (|2:29]) . 

Though we have thus demonstrated the feasibility of going to 4-loop level in thermal field 
theory, we would at the same time like to stress that the present computation was based on a 
"brute force" approach for the evaluation of the 3-loop and 4-loop sum-integrals (cf. Appendix 
A). Such an approach requires a lot of patience, and is simultaneously susceptible to errors. 
Thinking about the case of QCD, where a much larger set of genuine 4-loop integrals will 
appear, it would clearly be most desirable to develop somewhat more automatised techniques 
for the evaluation of the loop integrals, in analogy with what has been achieved in recent 
years for massive vacuum integrals at zero temperature [28\ [29], in order to carry out the 
computation in a more controllable fashion. 

Finally, we should point out that the model we have considered has also been a pop- 
ular testing ground for many different theoretical tools as well as improved approximation 
schemes (30j-|47j. in the latter case with the goal of learning something about their numerical 
convergence. We have shown here that, for N = 16, the resummed weak-coupling expansion 
does appear to show some convergence up to the point where the scale-invariant effective 
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coupling constant ^3, defined in Eq. (j3.15p . reaches values ~ 4.0. This includes all val- 
ues relevant for QCD above the temperature T ~ A^g [16], and therefore should perhaps 
be interpreted as a positive feature from the point of view of the applicability of resummed 
perturbation theory at high temperatures, even though it of course must be stressed that 
such a direct comparison between scalar field theory and QCD is far too naive. 
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Appendix A. Details of the computation 
A.l. Self-energies 

A.1.1. n(p) 

Before attacking the two sum-integrals in Eqs. (j2.2ip . (j2.22p . let us review a few straight- 
forwardly verifiable properties of the functions n(P), and n(P), defined in Eqs. (j2.8p . (j2.9p . 
[We will also define a third similar function in Sec. IA.1.21 denoted by n(P).] As has been 
shown in Ref. the first of these can be written in the form 

n(p) = n(o)(p) + -^ +An(p), (A.i) 

where 

n(o)fp) = ^ = mrHi-e) 1 

^ ' " (4^)2-e r(2-2e) (P2)e ^^'^^ 

denotes its zero-temperature limit and is responsible for its leading UV {P^ 00) behaviour. 
The leading UV behaviour of the rest, i.e. the finite-temperature part of the function, is on the 
other hand obtainable by letting P become arbitrarily large in either of the two propagators 
of n(P) and then integrating over the other, which produces 2Zi/P'^. Subsequently, the 
remaining part of the function, denoted here by An(P), behaves in the UV as 1/P^. 

At e = 0, the UV-finite subtracted function An(P) can be written in a simple form by 
using a three-dimensional (spatial) Fourier transform. This gives fT] 

An(P) = / d^rl e^P re-l-ol^ (coth r - J - J ) , (A.3) 
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where r = lirTr. 

Sometimes we need to refer to the finite-temperature part of n(P), defined as 'nP"\P) = 
n(P) — I1^'^\P). Its Fourier representation after setting e ^ is obtained from Eq. (jA.3j) 
by leaving out the last term — f/3 from inside the parentheses. In one instance, we will also 
require another version of An(P), in which its leading UV behaviour is further subtracted. 
This function is defined as 

A'n(p) ^ An(p) - ^ - (4 - 2^)^ ) (i - ^,o) , (a.4) 

with the Kronecker-(5 ((5pQ = (5pp^o) introduced in the last term for convenience. Its Fourier 
representation after e — > is equivalent to Eq. ()A.3P apart from having an extra term of the 
form + (1 — (5pQ)f'^/45 inside the parentheses. The constant Ji is evaluated in Ref. [1] and 
reads 

A.1.2. n(p) 

It will be convenient in the following to add a third function to those given in Eqs. (|2.8p . 
(j2.9p . We define this function through 

where II^'^\Q) denotes the zero-temperature part of n((5) (cf. Eq. ()A.2|1 ). 
An analogous reasoning as for n(P) produces for n(P) the representation 

n(p) = n(°)(p) +Jin(°)(p) + An(p) , (a.7) 

where 

and An(P) behaves in the UV like The only difference with respect to the previous 

calculation is that when taking the UV limit of the finite-temperature part of 11 (P), the 
two "propagators" of which 11 (P) is composed are not symmetric, and the dominant — > 
oo behaviour is obtained when the large momentum is routed solely through the function 
n(°)(Q), while the argument of the 1/(Q - Pf propagator is integrated over. 

As with n(P), we next derive a spatial Fourier representation for An(P). We will carry out 
the Fourier-transforms strictly in three dimensions, even though the functions transformed 
may contain e 7^ 0; this is sufficient since, as we will see, no divergences in e appear in the 
formal coordinate-space representations (there are in fact divergences hidden in the lower 
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ranges of the r-integration in individual terms, but they cancel for the finite quantity An(P) 
that we are ultimately interested in). Defining the (inverse) Fourier transform 



d^c 



2-l/2-a^--3/2 



r 



3/2-a 



(27r)3 (g2 + g2)a Y{a 

and making use of Eqs. (lA.2p . (jA.Sp . we obtain the Fourier-representations 
1 



K3/2~a(kok) , (Ai 



{Q-P? 



nW(Q) 



n(°)(g) 



ld're''-^-P>^J^{qo-po,r;l) 



(^3^ gi(p-q)-rg-|go-poh 



1 

47rr ' 



, r2(l-e)r(e) _ , 



(47r)2-^ r(2 - 2e) 



d^^re*'!'" 
d^re*'!'" 
d^re*'!'' 



(47rr)3 

A^^ r3(l-e)r(-l + 2e) 



(l + |(Zo|r) + 0(e)| , 

.F(go,r;-l + 2e) 



(A.IO) 



(A.ll) 



(47r)4-2^ r(3 - 3e) 
2 



^ ld\e^^-'^!^e-M^-j^{3 + 3\qo\r + qy)+0{e)'^ , (A.12) 

where the symbol "~" is a reminder of the fact that the integrals are not well-defined around 
the origin, and we note that the divergent factors r(e) and r(— 1 + 2e) have cancelled against 
the corresponding ones in Eq. ()A.9p . 

Inserting now Eqs. (|A.10p . (jA.lip into the definition of n(P) in Eq. (|A.6p . produces 

n(P) « ^ yd^re^P-jl j;(l + |<7o|r)e-l''°l^'e-l«°-*'°l^' + 0(e)| 



(47r)4 

T 
(4^ 



d3re*P-'"<!e-|P"l''i 



rcsch^r + (2 + |pok)(|Po| + cothr) 



(A.13) 



with pq = pq/2ttT. Subtracting Eq. (IA.12P as well as a result obtained from Eq. (|A.lip . 



Tin(°)(P) 



(4vr) 



d3re*pr^e-|P°l^'^^(l + |po|r) + 0(e)|> , 



(A.14) 



leads to the well-defined form 

T 

e=0 



An(p) 



(47r)^ 



j dV-^ e*P'''e~IP°l''|fcschV+ (2 + |^o|r)(|^o| +cothi 



- \{3 + ?,\po\r +plr'^) -^{l + \pQ\r) 

T o 



(A.15) 



where the integration is convergent around the origin such that we have replaced the symbol 
"~" with equality. We note that the expression inside the curly brackets in Eq. (jA.lSh 
behaves at small f as f^, and An(P) consequently indeed vanishes at zero temperature, and 
behaves as at large P. 
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A.1.3. n(p) 

Finally, we divide the function n(P) into two parts by separating from it the contribution of 
the Matsubara zero-mode, 

= n,(p) + no(P), (A. 16) 

where the prime in the upper right corner of the sum-integral symbol signifies the leaving 
out of the zero-mode go = from the corresponding Matsubara sum. Both of these parts are 
UV-finite. In three dimensions, one obtains the representations 

MP) = ^ /d'rie^P■r5]i^e-(l-W*'--l)^ (A.17) 

^ T f^l ( 2 1 \{p-q?+Pn 

no(P) = -77^ / dg-o 7^ + TT- In 



(27r)2 7o |^P2 ' 2pq 

T \m\ T 



47r P4 2(47r) 



{p + qY + pI 
y"d3re^P-'"e-|P«l'' , (A.18) 



where in the latter case we have first (in d = 2> — 2e dimensions) subtracted from the integrand 
the term the integral of which vanishes in dimensional regularization but which 

renders the expression for no(P) IR convergent. 

A. 2. Strategy for determining 5i 

As often in the evaluation of multi-loop (sum-)integrals, the first and in some sense also 
the most important step is to divide the integrand into two types of terms: ones that are 
divergent, but sufficiently simple to allow for an analytic evaluation, and others that are 
perhaps complicated but both IR and UV convergent. In the case of the sum-integral 5i 
defined in Eq. ()2.2ip . we do this by first decomposing n(P) according to Eq. (jA.ip . and then 
interchanging the order of sum- integrals, P ^ Q, in the terms involving n('^)(P) and 2Ii/P^. 
This leads to the following decomposition of the original sum-integral: 

5i = ^A^(p)^,(p) + ^A^(p)^o(P) + 2Jl^^^ +^^^ , 

= + + + . (A. 19) 

In 5}^ we repeat the split-up procedure by using Eq. ()A.7p for n(P). Furthermore, we take 
care of IR divergences by separating the contribution of the zero-mode from the Matsubara 
sums where necessary: 

Si' = ^\u{P)Uo{P) + TA^^I^-^U{po = 0,p)Uo{po = 0,p) 
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= + (A.20) 

= + 5™'^ + 5™'% (A.21) 




_ ^IV.a ^ ^IV,b ^ ^IV,c ^ ^IV,d ^ ^^ 22) 

The only further manipulations performed here are the dropping of terms that vanish in 
dimensional regularization; this happens in Matsubara zero-mode integrals without scales, 
after writing AH ^ H - n(o) - 2Ii/P^ in S^^'^, n(^) U - n(o) in 5™'^, and AH ^ 

n_n(o)-Tin(o) in 

We have now separated the original sum-integral into ten pieces, which fall into the follow- 
ing sub-categories: 

• Finite terms that can be evaluated numerically: Sl, S^'^, S^^'"^, S^'^. 

• Possibly divergent terms that can be evaluated analytically through the introduction 
of l^eynman parameters: b-^^ , b-^ , . 

• Terms that are trivial to compute due to Eq. (j2.19p : ^J^^''^, S^''^, S^''^. 

Below, we will go through the evaluation of the sum-integrals of the first category in detail 
and outline the calculation of those in the second category. 

A.2.1. Sl 

It is straightforward to see that S\ is finite in three dimensions, so we will immediately set 
e = 0. Expressing An(P) and Ilr{P) in terms of their Fourier representations, given in 
Eqs. (IA.3P and (lA.lTh . enables us to perform the p-integral in the definition of S\ to give a 
(5-function in coordinate space. This leaves us with the result 

7^3 f n 1 / _ 1 f\ sr^ sr^ f^-ilpol+llol+llo-poDr 



^ 2(47r)4 J \ f 3 J ^ ^ 



f°° 1 / 1 r\(l 
— — r / dr - ( cothr | < — cothr ln(l — e~^^' 

X 0.0269726622737(1) , (A.23) 



(47r)4 

where we have analytically performed the sums over po and qq (in this order) and later dropped 
the bars over the dimensionless coordinate variable. The last one-dimensional integral is of a 
form that might allow for an analytic evaluation, but for the purposes of this paper we have 
simply computed its value numerically. The number in parentheses estimates the uncertainty 
of the last digit. 
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A.2.2. 5"'^ 



For S^''^, we proceed along the lines of the previous section, employing this time the Fourier 
representation of JIo{P) given in Eq. (jA.lSh . After scaling all variables dimensionless, we 
obtain 

= -7^ /d^rl fcothr-i --We-2H'- 



ar coth r 



(47r)4 7o V r 3 y - 1 

X 0.0134942763002(1) , (A. 24) 



(4^)4 

where the last integration was performed numerically!! 
A.2.3. 5™'*' 

Expressing and n''"^-'(P) in terms of their respective Fourier integrals, we again perform 
the three-dimensional momentum integral and end up with one coordinate space integral and 
an infinite Matsubara sum. In dimensionless variables, we get 



g-2ln|r 



, .n 

dr I cothr 1 Infl — e~'^^^ 



3(4^)4 7o V 

rp2 



X 0.0625154109468(1) , (A.25) 



(4^)4 

where the last integration was performed numerically]! 
A.2.4. Sl^'"" 

To obtain Sj^''', we repeat the calculation of Sj^^^''*, but simply replace n(-^)(P) by An(P) 
and accordingly Eq. (jA.3p [without — r/3] by Eq. (jA.15p . This gives 

QlV,a _ — ^ f ^^3^ 'y^ — - i rcsch^r + (2 + |n|r)(|n| + coth r) — 
2(47r)^ / ^-^ \n\ 

1 r 1 I 

(3 + 3\n\r + n'^r^) - - (1 + \n\r) U^^l"!'' 

r 3 I 



^Its analytic value is T^[l + 7e + 71^/36 - ln(27r)]/(47r)*. 

^Its analytic value is T^(7e + 27i)/3(47r)*, where 71 refers to the first Stieltjes gamma constant (cf. the 
explanation following Eq. (|3.12p '). 



21 



dr -77 s -75 r coth r — 1 — — + 



(47r)4 Jq [e^^-l \ 3 

+^ ln(l - e"^'') ^3 + y - r^csch^r - 2r coth r 

~ ~ (47r)^ ^ 0.0004627085472(1) . (A.26) 
A.2.5. and 5^^'^ 

The sum-integrals iS"'*^, 5*"^'^ and S^'^ have one thing in common: the "outer" sum- integral 
is restricted to the Matsubara zero-mode po = 0, while the "inner" loop contains only one 
Matsubara sum. This implies that by introducing Feynman parameters, it is always possible 
to perform all of the momentum integrals analytically, leaving in the end a simple sum of 
a form that produces just a Riemann ("-function. As a first step, we quote the following 
results for the "static limits" (po = 0) of the various functions, obtained through Feynman 
par ametrization : 

nr„ -n„^ - t.\^^ '^'^'^^ r(3/2 + e)r(-i/2-6) i 

n = 0,;^ = /; ^ ' . > / dx x o 21 1/2+2. ' ^"^^ 

which we then substitute into the definitions of the sum-integrals in question. 
For SY'^, we obtain in the fashion described above 



,3 ,6e2^+'^V^ r(3/2 + 6)r(l/2 + e)r(-l/2 - e) 



,0 



(47r)3-2e r(-e) 
1 r d3-26p 1 



(27r)3-2e (p2)3/2+6[^(l _ ^)p2 + gg] 1/2+6 



^3.fe 2^+^^0Fr(3/2 + 6)r(l/2 + 6)r(-l/2-6; j ^,^2^-l/2-3. 
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where in the last form we have through rescalings of integration variables reduced the sum- 
integral into a product of sums and integrals that can each be trivially evaluated analytically. 
Performing this task, the final result for the function reads 



22 



To evaluate ^J^^'*^, we note that a calculation exactly parallel to the above produces 



III,b 



9T ^2.4. r(l/2 + 6) 



d3 



(47r)3/2- 
r 27r3/2-^ 1 



(2^)3-2. ^4[^(l _ ^)p2 + g2]l/2+e 



lr(3/2-e) (2^)3-2^ J, "^^(^2 + 1)1/2+ 



-2-2e 



(4^)4 36 
For S'J^''^, we finally get 

/3r2A4^ r(-l/2 + 2e) 



(A.32) 



^IV,b 



(4^)3/2- 



(l-x)- 



/3r2A4^ r(-l/2 + 2e) 
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A.2.6. 5™'% 5^^'^= and 



jIV,d 



The last three sum-integrals are trivial to evaluate, as their analytic values can be obtained 
by a straightforward application of Eq. (12.191) . The sum of the three gives 
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where 71 refers to the first Stieltjes gamma constant (cf. the explanation following Eq. (j3.12p ). 
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(A.34) 
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A.2.7. The full result for 

Collecting all the parts above, the final result for the sum-integral Si is seen to read 
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Evaluating the non-logarithmic terms numerically leads to the result in Eq. (|2.23|) . 



A. 3. Strategy for determining 52 

As with ^i, we begin the evaluation of our only genuinely 4- loop sum- integral 5*2, defined in 
Eq. (j2.22p . by dividing it into several pieces that are then calculated separately. Defining a 
regularized, or subtracted, zero-temperature limit for the function n(P) by 
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we write n(P) = l/(47r)^e + ns'^'*(P) + 11*^"^^ (P), and drop terms that vanish in dimensional 
regularization. This leads to 
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Here, we again decompose the various sum-integrals into further pieces according to 
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This time we classify the various parts as 
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• Finite (but comphcated) terms that can be evaluated numerically: S^^, S^'^^, S^^''^ . 

• Possibly divergent terms, whose divergent parts can be evaluated analytically and finite 
parts numerically: S2 ■ 

• Terms that are trivial to compute due to Eq. (j2.19p : 8^2'^ ■> S^2^'^ , 8^2^''^ 1 ^2^ 7 ^2- 
We move on to present the evaluation of the first four sum-integrals in detail. 

A.3.1. Sl'"" 

The term 82'^ is entirely finite in three dimensions, so we set e = and proceed to its 
numerical evaluation. There are two comparably simple ways of doing this, as one can either 
work entirely in coordinate space, or use the relation 



T 



n(^)(p) = 
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sm pr 



coth r — — 
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-\po\r 
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and perform the sum-integration over P directly in momentum space. Here, we choose the 
latter approach, but have verified the result also by the former method. Scaling all variables 
dimensionless, we straightforwardly get 
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(A.42) 



The summation can be accelerated by noting that for large n, the integral inside the curly 
brackets yields pj'iij?' +^1^), and the subsequent p-integral 7r/432n'^. This leading behaviour 
can be subtracted and the corresponding sum carried out analytically. 



A.3.2. 8. 



II,a 



The only potential divergence in ^2'^ is of an IR type, and produced by the zero-mode piece 
Po = of the 1/P^ subtraction term. However, there is in fact no divergence in dimensional 
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regularization, given that the zero-mode piece does not contribute: J d^~^'^p [n('^)(po = 0,p) — 
l/(47r)^e]/p^ = 0. We may therefore set e = from the beginning, if we agree to throw away 
the subtraction term for pQ = 0. Using furthermore the result of Eqs. (2.32) and (2.34) of 
Ref. flj, 
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as well as the representations for Ili'^^ and 11*^-^) as in Eqs. ()A.36p . ()A.4ip . we thus obtain 



II,a 



2 ^ roo 

- ^ y dp ln(p^ + n^) X 



3k /T_ 
(47r)4 1^12 



"'dr^ fcothr-ile-H^' 
r \ r J 



p- 



1 



16 



9 (p2 + n2)2 

C^(-3) 
C(-3) 



(A.44) 



For large n, the expression inside the curly brackets approaches Ap [p — 3n )/135(p + 
and the subsequent p-integral yields -7r[7 — 12 ln(2?7-)]/3240n^. This term can again be 



n 



2\4 



subtracted and the corresponding sums carried out analytically, to accelerate the convergence 
of the remaining numerical summation. We thus get 
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A.3.3. 52"'*" 

With ' we begin the calculation by noting that if we write 



(A.45) 
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and use the result in Eq. (D.14) of Ref. [1], 
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, (A.47) 
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as well as that in Eq. (|2J9]) of our Sec. H in order to show that ^pA'n(P) = ^p[n(P) - 
n(0)(p) -...] = 0(e), so that the K -term can be neglected, we can immediately reduce the 
evaluation of the original sum-integral into that of the simpler function 

In performing the latter task, we again have a choice between working entirely in coordinate 
space or combining (spatial) coordinate and momentum space techniques; this time we for 
illustration choose the former approach. 

We begin the calculation by recalling the result of Eq. (jA.9p . from which we obtain by 
differentiation with respect to a 



f d^p _ip.rr, , 2 , 2m2 e-H'- r ,^ , , ,/ ^ 2\m\ 

g|m|r 

— {1 — |m|r)Ei(2|m|r 



p— |m,|r 

f{r,m). (A.49) 



Here, we have used the relation 



in which the exponential integral Ei is related to the exponential integral function Ei through 
Ei(z) = — Ei(— z). We also note that for m = 0, the function f{r,m) becomes 

/(r,0) = 2(lnr-l + 7E). (A.51) 

Going now to coordinate space in the sum-integral and scaling the integration variables 
dimensionless, we obtain 

^P^'n(P)(ln^J = T^Epr^f(r,n)x 
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where all terms are both IR and UV convergent thanks to the (1 — (5„,) factor above. Given 
that it is perfectly finite, the zero- mode integral is furthermore simple to evaluate through 
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a procedure introduced in Ref. [T]: one temporarily introduces a convergence factor into 
the integrand, writes cothr = 1 + 2/(e^^ — 1), throws out all terms that are simple powers 
of r, performs the remaining integrals and finally proceeds to the limit a — > (which result 
can also be reproduced numerically). A straightforward calculation employing numerical 
integration for the non-zero mode terms then yields 



I^A'n(P) (in 
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Like before, the large-n behaviour can be worked out analytically and subtracted, to accelerate 
the convergence of the summation: it is 2T^[11 — 241n(2n)]/9072n^. 
Collecting everything together, we finally have 
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A.3.4. 5. 



i,b 



The function Sl^^ has one unique property: while entirely UV finite, it is the only sum-integral 
encountered in the present calculation that has a non-trivial logarithmic IR divergence that 
does not vanish in dimensional regularization. To see this, we note that the contribution of 
the Matsubara zero-mode to n(-^)(po = 0,p) [the superscript (T) is irrelevant here] reads 
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which leads to an IR singularity in the integral defining but cannot be separated from it 
without causing new UV divergences. 

To facilitate an analytic evaluation of the divergence, we add and subtract from the inte- 

grand of S2 a term of the form IIjj^ (po = 0,p) _ a'^^a 1 where the in principle arbitrary 
regularization mass m is chosen to be m = 27rT for computational convenience. This enables 
us to write 
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where the first piece has the virtue of being finite at d = 3, while the second one is straight- 
forward to evaluate analytically. Substituting here the Fourier representation of II^^^ (P) as 

(T) 

well as the above form of Iljj^ (po = 0,p), we obtain after scaling all variables dimensionless 
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The remaining five sum-integrals are again obtained through a straightforward application 
of Eq. ()2.19p . as they can be written in the forms 
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the sum of which gives 
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A. 3. 6. The full result for S2 

Taking the sum of all the various pieces of S2 displayed above, we can now finally write down 
the result for the entire sum-integral in the form 
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On the last row, we have separated the IR-singular contribution from 5*2' in Eq. ()A.57p . 
Evaluating the non-logarithmic terms numerically leads to the result in Eq. (j2.24p . 
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